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1. Introduction 

The Hurwitz zeta function, one of the fundamental transcendental 
functions, is traditionally defined (see by the series 

oo -, 

c(^'^) = E(^^' 3^(«)>o- (1) 

It admits an analytic continuation to the entire complex plane except 
for the simple pole s = 1. The Riemann zeta function (^(s) is a special 
case of z) 

C{sA) = C{s). 

The Hurwitz function has quite a few series and integral representations 
(see 15]). The most famous is the Hermite integral: 



, z^-' , ^ sin(garctan(f)) 
Cis,z) = — + 7 + 2/ — -dx, 



2 S-1 Jo (a;2 + 22)^/2(g2^x _ 1) 

s / l,5R(z) > 
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from which one can deduce many fundamental properties of the Hur- 
witz function, including the asymptotic expansion at infinity: 



l-s ^-s "1-1 



Z Z 



The Hurwitz function is closely related to the multiple gamma function 
Tn{z) defined as a generalization of the classical Euler gamma function 
r(2;), by the following recurrence- functional equation (for references 
and a short historical survey see [01 El ) : 

r„+i(z + l) = ^^±iM, zGC, nGJN, 

ri(z) = r{z), (3) 
r„(i) = 1. 

The multiple gamma function, originally introduced over 100 years 
ago, has significant applications in the connection with the Riemann 
Hypothesis. Montgomery |14j and Sarnak have conjectured that 
the limiting distribution of the non-trivial zeros of the Riemann zeta 
function is the same as that of the eigenphases of matrices in the CUE 
(the circular unitary ensemble). It has been shown in works by Mehta, 
Sarnak, Conrey, Keating, and Snaith that a closed representation for 
statistical averages over CUE oi N x N unitary matrices, when N ^ oo 
can be expressed in terms of the Barnes function G{z) = l/r2(z), 
defined by 



B2j r(2j + S-1) ^_2,-s+i , nr ^ 



G{z + l) = G{z)r{z), zeC, (4) 
G(l) = 1. 

Keating and Snaith conjectured the following relationship be- 

tween the moments of \C (^^ + itj \, averaged over t, and characteristic 
polynomials, averaged over the CUE: 

/I 1 /-^ 1 2A N 

log lim To / \Ci- + it)\ dt] = 

^\a{X)T^^ log^" {T)T Jo '^^2 ^' J 

A2(7 + l)-2Af: (-A)^2'=-l)|^^, 

k=2 '^^ ^ 

where a(A) is a known function of A G IN. The series on the right hand 
side of © is understood in a sense of analytic continuation (see P 
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for the method of evaluation of such sums), provided by the Barnes 
function 



|z| < 1, 

where 7 denotes the Euler-Mascheroni constant. Conversely, we find 
that the moments of |C + I in © are just a ratio of two Barnes 
functions: 

og lim To / \C{- + tt)\ dt]=-^^ ^— 

^ \a{X) T^'^ log^\T)T Jo '^2 ^' / G(2A + 1) 

The evidence in support of Keating and Snaith's conjecture is confirmed 
by a few particular cases A = 1,2,3,4, and the numerical experiment, 
conducted by Odlyzko [E], for T up to lO^^th zero of the Riemann 
zeta function. 

The identity © can be further generalized to 



and then evaluated in terms of the r-tiple gamma functions. For in- 
stance, with r = 3 we find 



E(--) 



k Cjk) 

+ 3z - 1 
12i 



4logr3(z + l) + ^logG(z + l)+ (7) 



+ 



7Z log(27r) 2C'(-1) 



\z\ < 1. 



In this paper we aim at developing a mathematical foundation for 
symbolic computation of special classes of infinite series and products. 
Generally speaking, all series (subject to convergence) of the form: 

00 

5^i?(i)log™P(j), mEiV, 
j=i 

where R{z) and P{z) are polynomials, can be expressed in a closed 
form by means of the multiple gamma function, which may further 
simplify to elementary functions. This algorithm also complements the 
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work previously started by Adamchik and Srivastava 1^ for symbolic 
summation of series involving the Riemann zeta function. The algo- 
rithm will expand the capabilities of existing software packages for 
non-hypergeometric summation. 



2. Asymptotics of C,'{—X,z) 

In this section, based on the integral ((21), we derive the asymptotic 
expansion of C'{t,z) = ^({t,z), when t = —A, A G Mq and z ^ oo. 
Differentiating both sides of (j2j with respect to s, we obtain 



C'{-X,z) = — —log 2-— log Z 



A + 1 ° 2 ° (A + 1 



,2 



^ ^^°o tan-Hf)cos(Atan-Hf)) ^^^ 
r°o log(x2 + 2;2)sin(Atan-i(f )) ^ 

Next, we expand the integrands into the Taylor series with respect to 
z. Taking into account that 

A/2 



^(Atan-(f))^E(-l)^(4)- 



and 

(x^ + z'f' sin (A tan-H^)) = f (-1)'= (^fc + 1^ 
we compute 

tan-^{^) cos(Atan-i(^)) 

^ ax = 



^A-2fc-1^2fc+l 



(e2^^-l)(x2 + z2)-A/2 



(9) 



and 
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log(x'' + z^) sin (Atan -"^(f)) 



dx 



A/2 
A:=0 



(10) 



(-1) 



2A; + 1, 



o27rx ]^ 



In the next step, we find asymptotic expansions of integrals on the 
right hand side of @ and (fTU]) when z — > oo. Expanding tan~^(|) and 
log (x^ + z^) into the Taylor series with respect to x and performing 
term by term integration, we arrive at 



x'" tan-i(^) 



1 



dx 



j=0 



4(2j + l)(j + fc + l) 



22r+3 



c27ra; i 



dx 



^ (-l)^.-^^-i?20-+.+l) 1 . 



4j(j + A: + 1) 



where Bk are the BernoulH numbers. Substituting these into © and 
(jlUp and the latters into ©, after some tedious algebra, we obtain 

PROPOSITION 1. The derivatives of the Hurwitz zeta function have 
the following asymptotic expansion when z ^ oo: 



^A+i ^A+l 

C'( — X,z) = log 2 logZ TT 

^ ^ A + 1 ^ 2 ^ (A + 1)^ 



+ E 



-2j+A+l 

"Mr 

z' 



2j-l 

^2, 5:(A: + 1) 

k=0 



2j-l 
/c + l 



(-A)^ 
1 



(11) 



-2j+A+lD 

logzv:^ ^^iL^hi^ + oi- 

^ fr{ (2j)! +^^^2r-A+l. 



where C'{t, z) = C{t, z), {-\)k = (-A)(-A + 1) . . . (-A + k-l) is the 
Pochhammer symbol, and [^] are the Stirling cycle numbers, defined 
recursively \lb^ by 



'n 


= {n-l) 


'n — Y 




'n — 


r 




'n 




k. 


k 


+ 


k- 












1. 




.0. 





1, n = 0, 
0, n / 



(12) 
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For A = and A = 1, formula (fTT|) yields the well-know asymptotics 

loe 1 
C'(0, N)=NlogN- -N + O(-), 

(13) 

/ iV^ 1 \ iV^ 1 1 

C'(-l,Ar)= !^ - tl + ±]iogN - — + — + 0(-^) 

In the similar way, we can derive the asymptotic expansions for higher 
order derivatives Q^'^^—X^N), where N ^ oo and A G 



3. Multiple gamma, zeta and the Hurwitz functions 

In [20], Vardi expressed the Tn function in terms of the multiple zeta 
function (n{s,z), as 

(-l)Nogr..(.) = -lim (M)-i: (-Iff : (14) 



k=l 



where 



and 



c„(., = E ■ ■ ■ E 



k=l 



OO OO -j^ 



.,=0 ,^^,{h + h + ... + K + zr 

The aim of this section is to find a closed form representation for 
logr„(2:) in finite terms of the Hurwitz zeta function, and vice versa. 

PROPOSITION 2. The multiple (z) may be ex- 

pressed by means of the derivatives of the Hurwitz zeta function: 



(_l)n n-L , s 

log r„(z) = ^ PkA^) ( C'(-fc) - C'i-k, z)] , > 0, 



k=0 



(17) 
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where the polynomials Pk,n{z) are defined by 

i-k-i (j - '^^ 



E (- 

j=k+i 



where 



are the Stirling cycle numbers. 



) 


'n' 




J. 



(18) 



The polynomials Pk,niz) can be envisaged as the generalized Stirling 
polynomials of the first kind, generated by 

71—1 re— 1 

Y[{k + X-Z) = Y,Pk,n{^)x^ 

k=l k=0 

with the following two alternative forms of representation 



Pk,n{^) 



and 



i=k+l 



kl dy^-^ (1 - y)i-^ 
z \(n-l)! 



.n — ij {i — 1)! 



k + l 



(19) 



(20) 



For z = 1 the polynomials Pk,n{z) are simplified to the Stirling num- 
bers: 



Pk,n{l) 



n — 1 
k 



(21) 



Interestingly, the polynomials Pk,n{z) were first considered by Mitri- 
novic il3j 40 years ago (with no relation to the multiple gamma func- 
tion) as a possible generalization of the Stirling numbers. 

The proof of Proposition 2 is based on the following two lemmas. 

LEMMA 1. The multiple zeta function C,n{s,z) defined by Mb]) may be 
expressed by means of the Hurwitz function 



n— 1 

^) = (r, - IV E ^J>(^) C(« 



■J,z) 



Proof. Recall the definition of the Stirling cycle numbers 

1 

Tv 

1=0 



k + n-V 
n-1 , 



in 



1 ! ^ 



(22) 
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Expanding /c* ^ = ((A; + z) — zf ^ by the binomial theorem, imphes 
that 



fc + n — 1 
, n - 1 



n i—\ 



in 



1 i-r 

1=0 j=0 \ / 



{k + zYi-zy-^-^ 



Interchanging the order of summation and making use of ()18|) , leads us 
to 



k + n-r 
n-1 , 



1 



(n-1) 



j=0 



Finally, we substitute this binomial coefficients representation into (|16|) : 



1 



(k + n-l 
^' n-1 



We complete the proof by evaluating the inner sum to the Hurwitz zeta 
function. □ 



LEMMA 2. The function Rn defined by il5\) may be expressed by 
means of the derivatives of the Riemann zeta function 



n-1 



Rn 



n 
k + 1 



("-l)'.=o 

Proof. Using Lemma 1 with z = 1 and formula l\'21\\. we get 



(23) 



(k-l) 



1)! ^ 



i=o 



k-l 



Cis-j) 



which, upon differentiation, yields 



dCk{s,l) ./f^^jr 

lim ^ = C 0) 5o,fc-i 

s^O OS 



{k 



k-l 

3 



C'(-j) (24) 



Here 5k^n is the Kronecker delta function. Next, we sum both parts of 
identity (|24l) with respect to k: 
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k=l 



k=j 



'(k-l)\ 



k-1 
3 



Observing, that the inner sum can be expressed via the Stirhng num- 
bers: 



E 



k=j 

we complete the proof. 



1 


-k-r 


1 


n 


{k-iy. 


. j . 


(n-1)! 


J + 1_ 



□ 



Proof of Proposition 2. The derivation of H17() immediately follows 
from formula dJ), by applying Lemma 1 and Lemma 2: 



iogr„(z) 



(n-1) 
1 



n-1 

n-j-l 



n - J 
k + l 



Interchanging the order of summation, we obtain 



iogr„(z) 



(-1) 



ra-l 



(n-1) 



WE -^.,n(^)C'(-J,^) + 



n—k 

c'i-k) y: (-1) 

i=o 



j=0 

z\ (n - 1)! 



J I {n-j- 1)! 



n-j 
k + l 



Finally, we note that the above inner sum is equal to representation 
(|20|) with a reversed order of summation. This completes the proof. □ 



In the second part of this section we invert Proposition 2 and express 
the derivatives of the Hurwitz function in finite terms of the multiple 
gamma function. The problem of representing derivatives in terms of 
other special functions was originally addressed by Ramanujan who 
studied functional and asymptotic properties of this function ^ 

cj)^{z)=C\-k,z + l)-C'{-k). 

It is unknown if Ramanujan was aware of Barnes' results regarding the 
multiple gamma function. 
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PROPOSITION 3. The derivatives of the Hurwitz zeta function may 
be expressed by means of the multiple gamma function Tn{z) 



C'(-n,z)-C'(-n) = (-ir^A:!Qfc,„(z)logrfc+i(z), > (25) 

k=0 

where the polynomials Qk,n{z) are defined by 

QkA^) = j:{l-zr''(''.){l] (26) 
j=k VJ / L 

and I are the Stirling subset numbers, defined by 

The immediate particular case of (|25j) is the formula for the triple 
gamma function: 

C'(-2,z)-C'(-2) = 21ogr3(z) + (3-2z)logG(z) (28) 

+ (1 -z)2logr(z), K(z) > 

For z = the polynomials Qk,n{z) simplify to the Stirling subset 
numbers: 



j=k 

LEMMA 3. Polynomials Pj^n{z) and Qk,j{z) satisfy the following dis- 
crete orthogonality relation: 

n—l 

E i-^y~''QkA^)PjA^) = h,n-l (29) 
j=k 

where 6k,n-i is the Kronecker delta. 

Proof. The proof is primarily based on the discrete orthogonality 
relation for the Stirling numbers: 

6. 



r ' 1 r ^ 
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□ 



Proof of Proposition 3. The proof readily follows from p5|) by re- 
placing \ogTk-\-iiz) by H17() and then using lemma 3. □ 

Polynomials Qk,n{z) generate a great deal of new integer sequences 
as well as closed form representations for the existing ones. Here are 
two examples from the online encyclopedia of integer sequences |18| . 

Example 1 

The sequence A021424: 1,16,170,1520,12411,96096 made of coef- 
ficients by X in the expansion of 1/((1 — x)(l — 3x)(l — 53;)(1 — 7x)). 
This sequence has a closed form representation 

2"-=^Q3,n(^) = ^(7" - 3 5" + 3"+i - 1) 

Example 2 

Riordan's sequence A000554 (labeled trees of diameter 3 with n 
nodes): 12,60,210,630,1736,4536,11430. This sequence is generated 
by the coefficients of in the polynomial 2 * Q2,n{^ — z), which is also 
can be written as 

52Q2,„(1 - Z) 











n - 2 ) ' 




2^0 




I 2 ; 



dz'^ 



4. Summation 

Consider the class of infinite sums 



J2Rik)logP{k) (30) 



k=l 



where R{z) and P{z) are polynomials. We will show that all such 
sums can be expressed in finite terms of the multiple gamma function. 
Clearly, (|30|) is a linear combination (subject to the branch cut of log) 
of 



N 



$7v(z) = Jim ^Flog(fc + z) (31) 
Thus, it is sufficient to find a closed form for H31j) . Expanding 
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F = {{k + x)- xf 
by virtue of the binomial theorem, we get 

T) / \ 



•^Niz) = lim ^ r i-iyzP-^F^iz 



3=0 ' 



where 



N 

FNiz) = J2ik + zy\ogik + z) (32) 

fc=0 

Assuming the analytic property of the Hurwitz zeta function, (|32|1 
can be written in the form 

, N , 

s-^oo ds ^ [k + zy^ 

and therefore 

p 



(33) 

(c'{-j,N + z + l)-C'{-j,z)] 



Further, we replace the first cumulant C'(~J) z) by a linear combination 
of the multiple gamma function (see Proposition 3), and the second 
cumulant (!^'{—j,N + z + 1) by asymptotic (lllj) . When the original 
series (|3()|) converges, the A^-dependent term will vanish as A ^ oo. 
However, if series (|3()j) diverges, formula (|33)) combined with (|llj) will 
provide a zeta regularization of (|3flj) in the Hadamard sense. 

In order to make this idea clear, we will consider a few examples. 
4.1. Dilcher's sum 

In [3| , Dilcher introduced a particular generalization of the Euler gamma 
function T(x) which is related to the Stieltjes constants 7^, as T{x) is 
related to the Euler 7 constant. The Stieltjes constants 7^ are defined 
as coefficients in the Laurent expansion of the Riemann zeta function 
C{s) at the simple pole s = 1: 

fc=0 
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The Euler 7 constant is a particular case of the Stieltjes constants 
7o = 7. Dilcher derived a few basic properties of the generahzed gamma 
function as well as its asymptotic expansion when z — > 00. The fol- 
lowing infinite series appears as a constant term in the asymptotic 
expansion: 

00 , , 

Dk = Y^ log^j + 2) - 21og'=i + log^j - 2) (34) 
i=i 

For k = 1 the sum (|34|1 simplifies to Di = log {2/tt). For other values 
of > 1, the closed form was unknown. In this section, we apply the 
technique of the multiple gamma function and evaluate the logarithmic 
series (|34|) in closed form. 

We begin with a finite sum 

N , , 

Dk = lim ^ log'^O' + -) - 2 log'^j + log*^(j - -) (35) 

and then rewrite it in terms of the Hurwitz function. The following 
chain of operations is valid (assuming the analytic property of the 
Hurwitz function) 

(-1)'= y iog\j + c) = lim ^ y ^ = 

1 ^ 1 \ _ (36) 

C\0,c+1)-C\0,c + N + 1) 
Applying this to (|35|) , yields 



Dk = -2CH0, 1) + C'(0, ^) + C'(0, ^) + 
lim (2 C^O, N + l)-CH0,N + l)-C''i0,N + l[ 



(37) 



Now we use the asymptotic expansion of C^i^^ N), N ^ 00. From the 
Hermite integral (j^J it follows that the dominant asymptotic term of 
C^{0,N) comes from the first two terms 



^ s ^1 s 



C{s,z) = — + 



2 s-1 

This leads to the following asymptotic expansion N ^ 00 
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J=o (38) 

where {k — j + l)j is the Pochhammer symbol. Thus, to prove that the 
limiting part in (|37|) vanishes to zero, it is sufficient to show that 



lim 



2(7V + l)logJ'(iV + l) 



iV + ^)log^(iV + ^ 



(39) 







for any j = 1,2,... ,k + l. This follows immediately from the fact that 

N ^ CO 



log{N + c) = logN+^ + 0^^ 

and the sum of coefficients by log''"^+^(iV) in ^ IS zero. 
Therefore, we have 



Dk = i-ir ( - log" 2 - 2 C'=(0) + 2 C"(0, -) 



(40) 



Note, C^(0, ^) in the right hand-side of (gOl) is a linear combination of 
^=(0): 

(^=(0, ^) = -ilog'^2 + J2 (^) C'-^(0) log^' 2 
j=i \3 / 

Combining this with 1)40^ . yields 

PROPOSITION 4. The infinite series ^341 ) can be expressed in finite 
terms of the derivatives of the Riemann zeta function 

^-^Dk = - log^ 2 - C'(0) + J2 C'-^(0) log^' 2 (41) 
In particular, for k = 2 and A; = 3 we have 



^2 



vr 
12 



+ log^vr - 3 log^2 -Y - 271 



(42) 
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= 3D2 log 2 + 2 C"'(0) + 9 log (2^) log2 2 + log^ 4 (43) 

respectively. 

Note that by the functional equation 

C{s) = 2^7r^-ir(l - s) sin (^) C(l - s), ^{s) + 1 

the Riemann zeta function admits a meromorphic continuation to the 
entire complex plane. Differentiating this with respect to s and setting 
s = 0, allows us to express derivatives C'^(O) i^i finite terms of the 
Stieltjes constants (see |2] for a closed form evaluation). Therefore, 
formulas like (|42j) and (|43j) can be used for numeric computation of 
C'^(O) as well as the Stieljes constants 7^. 

4.2. Melzak's Product 

Let us consider the infinite product 

and evaluate it in closed form. Originally, this product was computed 
by Melzak ^2] in a particular case when x = \ 

lim TT 1 + - = — (45) 

k=l 

Later on, P. Borwein and W. Dykshoorn [S] generalized ()45() and com- 
puted it in terms of Bendersky's function, defined by (see [26]) 

li22...n" = ri(n+l), 
which is closely related to the Barnes function G 



ri(n+l) 



n! 



G{n + 1) 



Note that 1)44(1 can be written in an alternative form by changing the 
product limit from 2N to 2N + 1: 

hm IT 1 + ^ = e^^ hm 17 1 + — 
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In this section, using the multiple gamma function technique, we reeval- 
uate product (|44() . derive its few particular cases in terms of known 
constants and then generalize it to 

2^ / 4x2 



1™ n 1 ,2 

PROPOSITION 5. The following identity holds for ^{x) > -i 



Proof. We first convert this product into a finite sum by applying 
the logarithm to it: 



JV— >CX3 ,\ k 
k=l 



(47) 

exp ( lim (-1)*'^ \ \ogk- log(/c + 2x) 



fc=l 



In the next step we split the finite sum in the right hand-side of (|47|) 
into three sums: 



2N 2N 



(-l)'''A;( log k - log{k + 2x)) = 2 ^ (-1)'' log(A; + 2x) + 

k=l k=l 
2N 2N 

J2 i-l)^klog k-J2 + 2x) log(A; + 2x) 



k=l k=l 



and evaluate each of them in terms of the Hurwitz function. For arbi- 
trary c such that 5R(c) > 0, we find 

2N r. 2N (_-,\k 

E (-i)'(t + iog(t + c) = - ta^ a; E ^ = 

iVl„g2 + 2C'(-l,^)-2C'(-l.£±l) + 
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and 

2N 



\og(k + c) = C (o. ^) - C (o. ^) + 

C'(0,^.A')-C'(0.^.A^) 

where derivatives C'(— A,2;) are understood by Now we use Propo- 
sition 1, in particular asymptotic expansions ((T^ . to get 



2N 

lim V(-l)''A;(logA:-log(A; + 2x 



k=l 

2xlogr(x + ^)-2xlogr(x + l)- ^^^^ 

2C'(-1,:e + ^) + 2C'(-1,x + 1)- 

log2 , 
^-^-3C(-1) 

In view of Proposition 3 with n = 1, we convert the derivatives of the 
Hurwitz function to the Barnes function: 

C'(-i, z) = - iogG(z + i) + 2 iogr(2) + c'(-i) 

Combining these with (|48|) . after some algebraic manipulations, leads 
to (jini). □ 

Here are a few particular cases of H46() : 

lim IT (1 + t) = (49) 

iv-.c»ii/ 2fc^ r(i) ^ ^ 

where A is the Glaisher-Kinkehn constant and G is Catalan's constant 
defined by 
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iog^ = Y^-C'(-i) 



fro (2^ + 1) 



respectively. 

The Melzak product can be further generahzed. Here is one of such 
formulas 

PROPOSITION 6. The following identity holds for 3?(x) > -i/ 



fe=l 



2x 



TT 



G(l + x)G(l -x) 



i r3(i-x)r3(i + x) 



(52) 



We skip the proof of this proposition, since it could be done in the same 
manner as in Proposition 5. 

CORROLARY 1. Identity i[5ljj) can he further simplified to 



hm rr 1 - — 



tan(7rx) exp 



2x\l + 7^i) + ^^ + 



(53) 



^ (^Li2(a;) - Li2(-u;)^ - (^h{(^) - ^h{-i^] 
where ui = exp(27ria;) and Lifc(u;) is the polylogarithm, defined by 



Lifc('^) = -r, fc > 1, |a;| < 1 



fc=o 



Proof. The proof follows straightforwardly from H52|) upon employing 
the reflexion formulas for the multiple gamma function: 



G{l+Z) I TT 

log — — = z log 



G{l-z) 



smvrz 



^-^B,{z) + ^U,{e'--) 
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2 log r3(l + z) r3(l -z))+ log G(l + z) G(l - z) 



z^logf-^) -T:izB2{z) + ^B^{z) + 

where Bj.{z) are the Bernoulli polynomials. □ 
Proposition 6 yields the following particular cases: 

,■ „ 1 /I 2G 7C(3)\ ,,,, 

k=l 

Formula follows from with x = \^ uj = i and taking into 
account 

Li2(i) -Li2(-i) = 2iG 



i TT^ 

Li3(i) -Li3(-i) = — 
Similarly, (fS5|) follows from with x = ^, u; = — 1 and 

Li2(l)-Li2(-1) = ^ 



Li3(l)-Li3(-1) = ^ 



tan ^^^(ttx) vr 



1™ 1 .9 
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